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PRELIMINARY  RESULTS  OF  AN  ANALYSIS  OF  BEAM 
NOISE  IN  THE  MEDITERRANEAN U) 


I.  INTRODUCTION 

(S)  The  beam  noise  measured  on  the  LAMBDA  II  array  system  during 
Task  V of  the  Mediterranean  ASW  Augmentation  Program  is  currently 
being  analyzed  in  two  30  Hz  frequency  bands:  a 275-305  Hz  band, 

measured  on  the  high  frequency  (HF)  array  segment,  and  a 30-60  Hz 
band,  measured  on  the  medium  frequency  (MF)  array  segment.  This 
report  presents  preliminary  results,  obtained  at  an  intermediate 
stage  of  the  analysis,  along  with  a description  of  the  metho- 
dology and  the  processing  algorithms. 

(C)  The  objective  of  the  analysis  is  to  statistically  characterize 
the  narrowband/narrowbeam  noise  over  broad  regions  of  the  frequency- 
azimuth  plane.  To  this  end,  the  total  beam  power  is  decomposed  into 
two  components:  a spectral  line  component,  which  consists  of  the 
narrowband  lines  which  are  resolved  in  frequency,  and  a background 
component,  which  consists  of  the  total  contribution  from  those 
sources  which  are  not  resolved  in  frequency  and  may  or  may  not  be 
resolved  in  azimuth.  The  spectral  line  component  determines  the 
false  target  environment  against  which  target  classification  algor- 
ithms must  operate.  This  component  will  be  characterized  in  terms 
of  the  probability  distributions  of  the  number  of  lines  in  specific 
regions  of  the  frequency-azimuth  plane  along  with  the  distributions 
of  the  spectral  line  power  in  these  regions.  The  background  com- 
ponent determines  the  minimum  level  that  a narrowband  line  must 
exceed  in  order  to  be  detected.  This  component,  together  with  the 
total  beam  noise,  will  be  characterized  in  terms  of  probability 
distributions  computed  as  functions  of  both  frequency  and  azimuth. 
The  results  of  this  analysis  will  provide  a statistical  basis  for 
evaluating  the  surveillance  capability  of  systems  with  frequency 
and  azimuth  resolution  characteristics  similar  to  the  LAMBDA  II, 

MF  and  HF  array  segments  in  high  shipping  density  regions  such  as 
the  Ionian  Basin. 

(S)  The  preliminary  results  in  this  report  are  based  on  the  time 
series  of  the  frequency-azimuth  distribution  of  the  total  beam 
power  and  its  spectral  line  and  background  components.  The  charac- 
teristics of  these  time  series  are  illustrated  in  Section  2 in  terms 
of  two  samples  from  both  the  30-60  Hz,  MF  band  and  the  275-305  Hz 
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(S)  HF  band.  In  Sections  3 and  4,  these  characteristics  are  interpreted 
in  the  light  of  a shipping  density  estimate  taken  during  the 
observation  period  and  their  implications  on  target  detection  and 
clutter  are  discussed.  The  processing  algorithms,  including  the 
procedure  used  to  reduce  the  effect  of  the  missing  hydrophones  on 
the  beam  pattern  sidelobe  structure,  are  discussed  in  the  Ap- 
pendices. 

(s)  The  preliminary  results  to  date  are  as  follows: 

The  beam  noise  in  the  HF  band,  as  well  as  the  MF  band,  is 
primarily  due  to  shipping. 

The  azimuthal  variation  in  the  beam  noise  in  both  bands  is 
in  qualitative  agreement  with  shipping  survey  data. 

Different  noise  sources  contribute  to  each  band,  indicating 
that  at  least  two  (acoustic)  classes  of  ships  are  present. 

The  observed  results  in  the  HF  band,  together  with  the 
results  of  ref  [1],  suggest  that  the  300  Hz  line  of  a 
Type  II  target  should  be  visible  in  an  essentially  clutter 
free  environment  approximately  50%  of  the  time 

The  50  Hz  line  of  a Type  II  target  should  rarely  be 
visible  in  the  MF  band  in  a low  clutter  environment. 

The  error  in  the  beam  noise  estimates  induced  by  the 
three  missing  hydrophones  in  the  HF  array  segment  can 
be  reduced  significantly  with  a simple  modification  to 
the  beamforming  algorithm. 

2.  THE  BEAM  NOISE  TIME  SERIES 

(C)  The  beam  noise  in  the  30-60  Hz,  MF  band  is  highly  anisotropic 
with  individual  samples  exhibiting  broadband  components  consisting 
of  both  the  narrowband  lines  and  the  frequency  continuous  contribu- 
tions from  the  noise  source  field.*  Furthermore,  the  extent  of  the 
anisotropy  varies  considerably  with  time.  During  approximately 
two-thirds  of  the  time  series,  a few  strong  noise  sources  dominate 
the  beam  noise,  causing  significantly  larger  azimuthal  variations 
than  during  the  rest  of  the  time  series  when  no  dominate  noise 
source  are  evident 


*In  this  report,  the  term  "broadband  component"  refers  to  a con- 
centration of  beam  power  across  the  frequency  band  at  a given 
azimuth.  The  term  "frequency-continuous"  refers  to  that  constituent 
that  is  smooth  with  respect  to  frequency. 
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(C)  The  beam  noise  in  the  275-305  Hz,  HF  band  is  similar  in 
character  to  that  in  MF  band  but  differs  significantly  in  the 
extent  of  the  azimuthal  variation,  the  number  of  narrowband 
lines  and  the  power  in  those  lines.  Furthermore,  there  is 
little  correlation  between  the  broadband  components  in  the 
HF  band  and  those  in  the  MF  band  indicating  that  the  noise 
sources  which  dominate  in  one  band  are  not  necessarily  dominant 
in  the  other. 

(U)  To  illustrate  these  characteristics,  two  samples  from  each 
frequency  band  have  been  selected  from  the  time  series.  The 
first  sample  is  chosen  from  a period  when  the  noise  in  the  MF 
band  is  dominated  by  two  broadband  components  while  the  noise 
in  the  HF  band  has  several  relatively  small  broadband  components. 
This  sample  is  fairly  typical  of  the  majority  of  both  the  MF  and 
HF  time  series.  The  second  sample  is  chosen  from  a period  when 
there  are  no  dominant  braodband  components  in  the  MF  band  but 
several  relatively  large  broadband  components  in  the  HF  band. 

This  sample,  although  less  typical  than  the  first,  illustrates 
the  extent  of  variation  in  the  time  series. 

2.1  The  Total  Beam  Power 

(S)  The  time  series  of  the  beam  noise  were  obtained  over  a 
forty-three  hour  period  during  which  the  LAMBDA  II  system  was 
towed  over  short,  east-west  tracks  in  the  southern  part  of  the 
Ionian  Basin.  During  this  period,  a source  transmitting  cw 
signals  at  280  and  40  Hz  with  source  levels  of  173  and  150 
dB//uPa,  was  towed  in  an  easterly  direction.  The  beam  power  was 
computed  using  a 10  minute  average  (21  time  samples),  at  .2  Hz 
intervals  in  both  the  MF  and  HF  frequency  bands.  The  broadside 
beamwidth  at  the  highest  frequency  in  each  band  is  approximately 
2.6  degrees;  this  value  increases  by  a factor  of  two  over  the  MF 
band,  but  only  by  a factor  of  l.l  over  the  HF  band.  Additional 
details  can  be  found  in  Appendix  A and  reference  [1]. 

(C)  The  first  samples  of  the  total  beam  power  are  illustrated  in 
Figures  2.1  and  2.2  in  an  intensity-modulated  format.  In  these 
plots,  eight  shades  are  used  to  represent  a 24  dB  range  in  4 dB 
increments  with  increasing  shade  darkness  corresponding  to  in- 
creasing beam  power.  The  black  shade  represents  beam  powers  in 
excess  of  the  78  dB  maximum  and  the  white  shade  represents  values 
below  the  54  dB  minimum.  All  values  are  relative  to  the  intensity 
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(C)of  a luPa  planewave  measured  at  the  hydrophone  inputs.  The  frequency 
axis,  which  spans  the  30  Hz  range  for  each  band,  has  tic  marks 
separated  by  1.76  Hz.  The  azimuth  axis  is  linear  in  the  sine  of, the 
azimuth  variable  with  zero  degrees  corresponding  to  due  north. 

(C)  The  total  beam  power  in  the  MF  band,  Figure  2.1,  is  dominated 
by  two  broadband  components  at  approximately  -30  deg.  and  -3  deg. 

These  components  yield  variations  with  azimuth  of  up  to  34  dB 
depending  on  frequency.  An  examination  of  these  components  at  the 
higher  frequencies  (narrower  beamwidths) , suggests  that  each 
component  is  due  to  more  than  one  noise  source.  This  is  born  out  in 
the  time  series  in  the  vicinity  of  the  sample  which  shows  that  each 
component  is  primarily  due  to  two  sources  moving  in  opposite  direc- 
tions. Also  evident  in  Figure  2.1  is  a harmonic  component  at  -75 
degrees  with  lines  spaced  approximately  6.5  Hz  apart  and  a smaller 
component  at  40  degrees  which  has  a peak  in  its  spectrum  at  approxi- 
mately 46  Hz.  The  weak  component  at  90  degrees  may  be  due  to  noise 
radiated  from  the  array  tow  ship. 

(C)  The  plot  of  the  beam  noise  in  the  HF  band,  Figure  2.2,  shows 
a very  strong  broadband  component,  several  lesser  components  and  two 
weak  components  in  addition  to  the  280  Hz  signal  line.  The  strong 
component  at  70  degrees  and  the  intermediate  component  at  10  degrees 
can  be  attributed  to  array  tow  ship  noise,  with  the  strong  component  due 
to  the  direct  path  arrival  and  the  intermediate  component  due  to 
arrivals  which  have  been  reflected  from  the  bottom  and  the  surface. 

This  hypothesis  is  support  by  the  fact  that  both  these  components  are 
present  throughout  the  time  series  at  the  same  azimuths  relative  to 
forward  endfire,  regardless  of  the  tow  ship  heading.  Furthermore, 
the  arrival  directions  associated  with  these  azimuths  are  consistent 
with  the  location  of  the  array  relative  to  the  array  tow  ship  and  the 
bottom. 

(U)  The  remaining  broadband  components,  which  exhibit  azimuthal 
variations  of  the  order  of  10  dB  or  less,  can  be  attributed  to  sources 


The  beam  values  are  obtained  by  subtracting  the  theoretical  signal 
gain  (30. dB)  from  the  observed  values.  Alternatively,  the  beam  power 
can  be  referred  to  an  uncorrelated  noise  field  by  adding  25.95  dB 
which  represents  the  sum  of  the  theoretical  array  gain  (18.75  dB) 

U^d  a bandwidth  correction  factor  of  5.2  dB. 

The  array  headings  for  the  first  and  second  samples  were  90  degrees 
and  89  degrees  respectively. 
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(U)  other  than  the  array  tow  ship.  The  three  components  at  -30,  ~5  and  38 
degrees  have  counterparts  in  the  MF  noise  of  Figure  2.1  which  was 
measured  20  minutes  later,  and  thus,  could  possibly  be  due  to  the 
same  noise  sources.  The  remaining  components,  however,  have  no 
counterparts  in  the  MF  band. 

(U)  Finally  it  is  noted  that  there  are  two  low  noise  directions 
in  the  vicinity  of  50  degrees  on  either  side  of  zero  azimuth. 

This  is  typical  of  much  of  the  time  series  for  both  the  MF  and  HF 
bands  and  is  consistent  with  the  shipping  survey  data  which  shows 
that  most  of  the  ships  are  concentrated  in  a 100  degree  sector 
about  zero  degrees. 

(U)  The  second  sample  from  the  MF  band,  measured  approximately  20 
hours  after  the  first,  is  plotted  in  Figure  2.3.  In  contrast  to 
the  first  sample,  there  are  only  a few  instances  where  the  beam 
noise  exceeds  the  78  <1 13  upper  limit  of  the  plot  and  the  azimuthal 
variations  are  typically  of  the  order  of  12  to  16  dB.  The  broad- 
band components,  which  span  much  larger  azimuthal,  sectors  than 
those  in  sample  one,  can  be  attributed  to  the  contributions  from 
many  sources,  none  of  which  dominate.  The  comparatively  low  levels 
in  this  sample  are  typical  of  the  latter  one-third  of  the  time  series.* 

(C)  The  second  sample  from  the  HF  band,  measured  12  minutes  after 
the  second  MF  sample,  is  plotted  in  Figure  2.4.  In  this  sample, 
there  are  five  broadband  components  in  addition  to  the  two  com- 
ponents due  to  the  array  tow  ship  noise  and  the  280  Hz  signal  line. 

In  contrast  to  the  first  sample,  these  components  are  not  restricted 
to  50  degrees  on  either  side  of  zero  azimuth  and  they  show  a somewhat 
larger  azimuthal  variation.  Furthermore,  three  of  these  components, 
at  -70,  25  and  45  degrees,  occur  at  azimuths  where  there  are  no 
counterparts  in  the  second  MF  sample.  Finally,  it  is  noted  that 
the  beam  levels  due  to  the  array  tow  ship  noise  are  generally  higher 
than  those  in  the  first  sample.  This  increase  is  particularly  evident 
in  the  reflection  of  a 285  Hz  spectral  line  appearing  at  10  degrees 
azimuth . 

2.2  The  Spectral  Line  and  the  Background  Components 

(U)  The  background  component  represents  an  estimate  of  the  frequency 
continuous  portion  of  the  total  beam  power.  The  spectral  line  com- 
ponent consists  of  those  narrowband  lines  which  exceed  the  background 
by  at  least  2.6  dB.  The  two  components  are  computed  so  that  they  sum 
in  a linear  scale  to  yield  the  total  beam  power.  As  a result,  each 
component  can  be  interpreted  as  the  total  power,  given  that  the  other 

*These  comparatively  low  levels  might  be  due  to  reduced  shipping 
activity,  since  the  latter  one-third  of  the  time  series  occurred  on 
a Sunday  morning. 
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42-6G  dB//uPa;  17  Nov.  74,  0652X  (II) 
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(U)  component  were  zero.  The  threshold  constant  of  2.6  dB  was  chosen  to 
limit  the  probability  of  detecting  a line  when  none  is  present  to 
.004.  The  algorithm  used  to  decomposed  the  total  beam  power 
into  its  spectral  line  and  background  components  is  described  in 
detail  in  Appendix  B. 

(U)  The  spectral  line  component  for  the  first  sample  in  the  MF  band 
is  plotted  in  Figure  2.5.  The  format  is  similar  to  that  of  the  pre- 
ceding section  except  that  the  dynamic  range  spans  only  18  dB  in  3 dB 
steps,  and  the  frequency  range  is  reduced  by  3.1  Hz  due  to  end  effects 
in  the  decomposition  algorithm.  The  spectral  lines  associated  with 
the  broadband  components  identified  in  the  total  beam  power,  are 
clearly  evident.  The  harmonic  structure  associated  with  these  lines 
is  consistent  with  a fundamental  frequency  corresponding  to  either  the 
blade  rate  or  the  shaft  rate  of  merchant  ships.  In  addition,  other 
lines  appear  in  Figure  2.5  which  do  not  appear  in  Figure  2.1  due  to 
the  4 dB  step  size.  The  power  in  the  spectral  lines  ranges  from  56  dB 
up  to  as  much  as  35  iB  for  some  of  the  lines  associated  with  the  domi- 
nant broadband  components  in  the  total  power. 

(U)  The  spectral  line  component  for  the  second  sample  in  the  MF  band 
is  plotted  in  Figure  2.6.  The  power  in  the  specLial  lines  for  this 
sample  is  significantly  less  than  for  sample  one  since  there  are  no 
dominant  broadband  components  in  the  total  beam  power.  On  the  other 
hand,  a detailed  examination  of  Figures  2.5  ana  2.6,  shows  that  the 
number  of  distinct  spectral  lines  in  both  samples  is  approximately 
the  same;  258  in  the  first  sample  and  260  in  the  second.  These  results 
suggest  that,  for  threshold  constants  of  the  order  of  2.6  dB , the 
number  of  detected  spectral  lines  should  be  relatively  independent  of 
the  distribution  of  the  spectral  line  power. 

(U)  The  background  components  for  both  samples  in  the  MF  band  are 
plotted  in  Figures  2.7  and  2,8.  The  f requency-continous  spectra 
associated  with  the  broadband  components  in  the  total  power  are 
evident  in  both  the  background  components.  The  two  large  broadband 
components  which  dominate  the  total  power  in  the  first  sample,  yield 
azimuthal  variations  of  up  to  25  dB  in  the  corresponding  background 
component  of  Figure  2.7.  In  contrast,  the  azimuthal  variation  in  the 
background  component  of  Figure  2.8  seldom  exceeds  15  dB.  Furthermore, 
at  the  azimuths  where  the  spectrum  level  is  large,  there  is  more 
variation  with  frequency  evident  in  Figure  2.7  than  in  Figure  2.8. 

This  can  be  attributed  to  the  fact  that  the  background  component  for 
sample  two,  to  a greater  extent  than  sample  one,  represents  the 
contributions  from  many  noise  sources  so  that  the  frequency  character- 
istics of  individual  sources  are  less  apparent. 

(U)  The  plots  of  the  spectral  line  component  and  the  background 
component  for  sample  one  in  the  HF  band  are  shown  in  Figures  2.9  and 
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Fig.  2.5  (C)  — Spectral  line  component,  MF  band,  sample  one 
56-7/}  HR// uPa;  16  Nov.  74,  1110Z  (U) 
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Fig.  2.7  (C)  Background  component,  MF  band,  sample  one 
56-7-1  dl!//u Pa;  16  Nov.  74,  lllOZ(U) 
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(U)  2.10  respectively.  The  narrowband  lines  associated  with  the  direct 
and  the  reflected  arrivals  of  the  array  tow  ship  noise  appear  at  70 
and  10  degrees  in  Figure  2.9  and  the  frequency-continuous  contributions 
appear  at  the  same  azimuths  in  Figure  2.10.  A comparison  of  Figure  2.9 
with  Figures  2.4  and  2.5  shows  that  there  are  fewer  spectral  lines  in 
the  HF  band  than  in  the  MF  band.  In  particular,  excluding  the  lines 
associated  with  the  array  tow  ship  noise  and  the  280  Hz  signal  line, 
there  are  approximately  175  lines  in  the  HF  band  as  compared  with 
approximately  260  linas  in  each  sample  from  the  MF  band.  Furthermore, 
the  spectral  line  power  in  the  HF  band  shows  less  variation  than  in  the 
MF  band,  with  all  but  a few  lines  having  levels  restricted  to  a 9 dB 
range  between  42  dB  and  51  dB . 

(U)  The  background  component,  Figure  2.10,  shows  frequency  and 
azimuthal  variations  that  are  slightly  less  than  sample  two  in  the 
MF  band  and  significantly  less  than  sample  one.  Excluding  the  com- 
ponent due  to  the  direct  arrrival  of  the  array  tow  ship  noise,  the 
largest  azimuthal  contrast,  of  about  15  dB,  occurs  between  the  low 
noise  direction  at  -50  degrees  and  the  broadband  component  at  -30 
degrees.  For  azimuths  within  + 50  degrees,  which  excludes  the  low 
noise  direction,  the  beam  levels  vary  by  less  than  10  dB  at.  all 
frequencies . 

(U)  It  is  important  to  note  that  at  no  azimuths  can  the  beam  level 
be  attributed  to  wind  generated  noise  alone.  For  example,  if  the  wind 
noise  is  assumed  to  be  spatially  and  temporally  uncorrelated,  then  an 
omni-directional  spectrum  level  of  55  dB//uPa,  which  would  correspond 
to  the  observed  sea  state  one,  would  result  in  a beam  noise  of  33.5  dB. 
This  level  is  6.5  dB  lower  than  the  lowest  beam  level  observed  in  the 
background  component  of  Figure  2.10, 

3.  IMPLICATIONS  FOR  TARGET  DETECTIONS  AND  CLUTTER 

(U)  A narrowband/narrowbeam  surveillance  system  operates  on  the  beam 
noise,  to  detect  the  presence  of  signal  lines  and  to  classify  fhes° 
detections  into  false  target  lines  and  real  target  lines.  Signal  lines 
due  to  both  false  targets  and  real  targets  are  detected  when  they 
exceed  the  background  by  some  threshold  constant.  An  increase  in  the 
threshold  constant  reduces  the  number  of  false  target  lines,  but  it 
also  reduces  the  percentage  of  time  that  the  real  target  lines  are 


*As  was  the  case  in  the  MF  band,  the  number  of  spectral  lines  in 
both  samples  in  the  HF  band  is  approximately  the  same;  172  in  the 
first  sample  and  178  in  the  second. 
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Fig.  2.9  (C)  — Spectral  line  component,  HF  band,  sample  one 
42-60  dB//uPa;  16  Nov.  74,  1050Z  (U) 
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(IJ)  visible.  Thus,  the  choice  of  a threshold  constant  represents  a 

compromise  between  the  extent  of  the  c lut ter , environment , and  hence, 
the  magnitude  of  the  classification  problem,  and  the  probability 
of  detecting  the  real  target  lines.  In  this  section,  the  results  of 
Section  2,  together  with  those  of  reference  [1],  are  used  to  inves- 
tigate this  compromise. 

(U)  The  dependence  of  the  number  of  detected  lines  on  the  value  of  the 
threshold  constant  can  be  investigated  in  terms  of  the  ratio  of  the 
total  beam  power  to  the  background  power.  Plots  of  these  ratios 
differ  from  the  plots  of  the  spectral  line  component  in  the  preceding 
section  in  that  they  sum  in  a decibel  scale  with  the  background  com- 
ponent to  yield  the  total  power.  In  the  context  of  "signal"  line 
detections,  these  ratios  represent  a "s ign_g].-plus-background  to 
background  ratio"  which  is  denoted  here  by"  SB/B.  The  number  of 
spectral  lines  that  are  detected  using  a given  threshold  constant  c, 
is  obtained  from  the  plot  cf  the  SB/'B  component  by  a count  of  the 
number  of  times  SB/B  exceeds  c. 

(U)  The  SB/B  components  for  both  samples  in  the  MF  band  and  the  first 
sample  in  the  HF  band  are  illustrated  in  Figures  3.1,  3.2,  and  3.3. 

The  range  of  these  plots  spans  0 to  12  dB  in  2 dB  steps  and  values 
less  than  2.B  dB  have  been  deleted  so  that,  the  only  spectral  lines 
that  appear  are  those  that  appear  as  spectral  lines  in  Section  2. 

(U)  A comparison  of  the  SB/B  components  for  the  two  samples  in  the 
MF  band,  Figures  3.1  and  3.2,  shows  that  even  though  the  power  in  the 
spectral  line  component  for  sample  one  is  significantly  larger  than 
for  sample  two,  the  SB/B  ratios  are  approximately  the  same.  This  can 
be  attributed  to  the  fact  that  an  increase  in  the  spectral  line  power 
is  accompanied  by  an  increase  in  the  background  component  due  to  an 
increase  in  the  frequency-continuous  portion  of  the  spectrum  from  the 
dominant  sources. 

(U)  In  the  plot  of  the  SB/B  component  for  the  HF  band,  Figure  3.3,  it 
is  seen  that  the  3B/8  ratios  are  distributed  over  a smaller  range  than 
those  in  the  MF  band.  As  a consequence,  small  increases  in  the 


The  classification  of  the  spectral  lines  depends  on  the  difference 
in  the  line  properties,  center  frequency,  bandwidth,  stability  etc,, 
between  the  false  target  lines  and  the  real  target  lines. 

The  SB/B  ratio  is  related  to  the  signal-to-background  ratio,  S/B,  by 

SB/B  = 10  log  (1  + S/B) 
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Fig.  3.2  (C)  — SB/R  component,  MF  band,  sample  two 
0-12  dB;  17  Nov.  74,  0640Z  (U) 
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(U)  threshold  constant;  will  have  a greater  effect  on  the  number  of  detected 
lines  i.n  the  HF  band.  To  illustrate  this  effect,  we  have  computed 
the  average  number  of  detected  lines  over  both  samples  in  each  fre- 
quency band  for  different  values  of  the  threshold  constant.  The 
results  are  plotted  in  Figure  3.4. 

(U)  Figure  3.4  shows  a dramatic  difference  in  the  dependence  of  the 
number  of  detected  lines  on  the  value  of  the  threshold  constant  for 
the  two  frequency  bands.  For  example,  in  the  HF  band,  an  increase  in 
the  threshold  constant  from  2.6  to  4 dB,  reduces  the  number  of  lines 
by  90  percent,  and  only  one  line  is  detected  for  a threshold  constant 
of  8 dB.  On  the  other,  in  the  MF  band,  the  increase  from  2.5  to  4 dB, 
reduces  t.he  number  of  lines  by  only  53  percent  and  18  lines  are  de- 
tected with  a threshold  constant  of  8 dB. 

(S)  Figure  3,4  illustrates  how  increases  in  the  threshold  constant 
reduce  the  clutter  due  to  the  false  target  lines.  The  compromi.se 
between  the  extent  of  the  clutter  due  to  false  targets  and  the  per- 
centage of  time  that  the  real  target  lines  are  visible  can  be  seen 
by  considering  two  examples:  the  detection  of  a 300  Hz  line  with  a 

mean  gpurce  level  of  150  dB  and  a 50  Hz  line  with  a mean  source 
level'**  of  153  dB, 

(S)  Consider  the  300  Hz  line  first.  Figure  4.1,  reference  [1],  shows 
the  probability  distribution  of  the  280  Hz,  s ignal-to-noise  ratio 
observed  during  the  experiment.  From  this  Figure,  and  the  discussion 
of  Section  5.0,  ref  [I],  it  can  be  concluded  that  the  50  percentile  of 
the  SB/B  ratio  for  the  300  Hz,  target  line  is  approximately  4.75  dB.*** 
It  follows  from  Figure  3.4,  that  with  a threshold  constant  of  4 dB, 
the  300  Hz  line  will  appear  as  a spectral  line  in  a low  clutter 
environment  approximately  50%  of  the  time.  On  the  other  hand,  it  can 
be  seen  from  Figure  4.1  (ref  [1])  and  Figure  3.4  that  decreasing  the 
threshold  constant  to  2.5  dB,  will  greatly  increase  the  clutter  'with- 
out significantly  increasing  the  percentage  of  time  that,  the  300  Hz  line 
appears. 


'it 

The  number  of  lines  in  the.  HF  band  does  not  include  the  280  Hz  line 
nor  the  lines  that  occur  between  5 and  15  degrees  and  70  and  90  degrees 
since  at  least  some  of  these  lines  can  be  attributed  to  the  array  tow 
ship  noise. 

"it  if 

_tiThese  source  levels  are  typical  of  a Type  II  target 

"The  value  of  4.75  dB  was  obtained  by  equating  s/N  with  S/B, 
subtracting  23  dB  to  account  for  the  difference  in  source  levels  and 
adding  .97  dB  to  account  for  the  difference  in  the  analysis  bandwidth. 
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NUMBER  OF  SPECTRAL  LINES 


MF  BAND 


HF  BAND 


THRESHOLD  CONSTANT  (dB) 

Fig.  3.4  (C)  — Number  of  spectral  lines  vs.  threshold  constant  (U) 
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(S)  The  same  analysis  cannot  be  applied  to  the  50  Hz  line,  since  the 
signal-to-noise  ratio  was  not  measured  during  TASK  V due  to  the  low 
source  level  at  40  Hz.  It  was  noted  in  reference  [1],  however,  and 
confirmed  in  this  anaLysis,  that  the  40  Hz  signal  line  is  only  rarely 
visible  above  the  background  in  the  MF  band.  Since  the  40  Hz  source 
level  is  only  several  dB  lower  than  that  expected  from  the  50  Hz  line, 
we  may  conclude  from  Table  3.1  that  detections  of  the  50  Hz  line  must 
occur  in  a high  clutter  environment. 

4.  INTERPRETATION  OF  THE  BEAM  NOISE  IN  TERMS  OF  THE  SHIPPING 
DISTRIBUTION 

(U)  The  azimuthal  variation  observed  in  the  beam  noi.se  is  i.n 
qualitative  agreement  with  a sample  of  the  shipping  distribution  taken 
during  the  43  hour  observation  period.  This  sample,  Figure  2 of  ref- 
erence [2],  shows  211  ships  in  the  Ionian  Basin  during  the  period  when 
the  first  samples  from  the  MF  and  the  HF  frequency  bands  were  selected. 
Of  these  211  ships,  aLl  but  49  are  located  within  a 100  degree  sector 
centered  at  zero  degrees  and  all  but  four  of  the  remaining  49  are 
located  in  the  image  of  that  sector  centered  at  180  degrees.  Thus, 
the  comparatively  low  beam  noise  values  observed  outside  of  the  100 
degrees  sector  in  both  frequency  bands  are  not  surprising. 

(U)  The  shipping  distribution  within  the  100  degree  sector  includes  a 
rather  broad  shipping  lane  that  traverses  the  Ionian  Basin  approximately 
250  miles  north  of  the  array  tow  site,  plus  a concentration  of  ships 
along  the  southeastern  coasts  of  Sicily  and  Italy,  300  miles  to  the 
northeast.  These  ships  are  sufficiently  dense  in  azimuth  so  that 
several  appear  on  each  beam  within  the  100  degree  sector.  It  is  the 
noise  from  these  ships,  both  narrowband  and  frequency-continuous,  that 
probably  accounts  lor  much  of  the  noise  in  the  background  component. 

(U)  In  addition  to  these  concentrations,  the  shipping  sample  shows  a 
scattering  of  ships  which  place  several  within  75  miles  of  the  array 
tow  site.  It  is  these  nearby  ships  which  may  be  the  sources  of  the 
dominant  broadband  components  which  appear  i.n  the  beam  noise  in  the 
MF  band  and  to  a lesser  extent  in  the  HF  band. 

(C)  The  interpretation  in  the  preceding  paragraphs  is  necessarily 
qualitative  since  the  shipping  survey  did  not  provide  any  information 
from  which  to  infer  the  source  level  characteristics  of  the  ships. 

The  need  for  this  information  was  seen  in  the  samples  of  the  beam  noise 
which  showed  that  many  of  the  broadband  components  that  are  signifi- 
cant in  one  frequency  band  do  not  appear  in  the  other  frequency  band. 
This  result  suggests  that  the  ships  must  be  divided  into  at  least  two 
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(C)  classes,  those  which  radiate  noi.se  primarily  at  the  higher  frequencies, 
such  as  fishing  vessels,  and  those  which  radiate  primarily  at  the  lower 
frequencies,  such  as  merchant  ships.  Thus,  without  knowing  the  spatial 
distribution  of  each  class  of  ships,  the  prediction  of  beam  noise 
statistics  using  a shipping  noise  model  must  be  viewed  as  highly 
speculative. 

5.  SUMMARY 

(U)  The  beam  noise  in  the  MF  band  is  highly  anisotropic  with  individ- 
ual samples  exhibiting  broadband  components  consi.r '-.ing  of  both  the 
narrowband  lines  and  the  frequency-continuous  contr i.butions  from  the 
noise  source  field.  The  power  i.n  both  the  spectral  lines  and  the  back- 
ground component  varies  considerably  over  the  43  hour  observation 
period.  However,  the  number  of  spectral  lines  and  the  signal-to- 
background  ratios  for  those  lines  are  relatively  stable. 

( 3)  The  beam  noise  in  the  HF  band  is  similar  in  character  to  that  in 
the  MF  band,  differing  primarily  in  the  extent  of  the  azimuthal 
variation,  the  number  of  spectral  lines,  the  power  in  those  lines, 
and  the  relative  contribution  of  the  noise  generated  by  the  array 
tow  ship.  The  number  of  spectral  lines  in  the  HF  band  is  approxi- 
mately 70  percent  of  that  in  the  MF  band  with  a threshold  constant 
of  2.6  dB.  However,  both  the  power  and  the  si.gnal-to-background 
ratios  of  these  lines  are  distributed  over  a much  smaller  range. 

As  a consequence,  the  extent  of  the  line  clutter  due  to  false  targets 
is  a much  more  sensitive  function  of  the  threshold  constant  in  the 
HF  band.  This  result,  together  with  the  probability  distribution  of 
the  signal-to-noi.se  ratio  from  reference  [1],  suggests  that  the  300 
Hz  line  of  a TYPE  II  target  will  be  visible  in  a low  clutter  environ- 
ment approximately  50  percent  of  the  time,  but  that  the  50  Hz  line 
will  almost  never  be  visible  in  other  than  a high  clutter  enviroment. 

(C)  The  beam  noise  in  the  HF  band,  as  well  as  the  MF  band,  is 
primarily  due  to  shipping,  even  in  the  low  noise  directions,  as 
evidenced  by  a comparison  of  the  observed  beam  noise  levels  with  that 
expected  from  wind  generated  noi.se.  The  azimuthal  variation  i.n  the 
beam  noise  is  in  qualitative  agreement  wi.th  that  expected  from  con- 
current ship  location  data  which  indicates  that  almost  all  of  the 
ships  are  concentrated  i.n  a 100  degree  sector  about  due  north  of  the 
nominal  array  location  and  the  image  of  that  sector  about  due  south. 

It  was  oDserved,  however,  that  there  are  many  instances  where  different 
noise  sources  contribute  in  each  band,  indicating  the  presence  of 
more  than  one  (acoustic)  ship  class.  As  a consequence,  the  knowledge 
of  the  spatial  distributions  of  each  of  the  ship  classes  may  be  re- 
quired for  an  accurate  prediction  of  the  beam  noi.se  statistics  using 
shipping  noi.se  models. 
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(U)  Finally)  i.n  Appendix  C5  it  is  shown  that  if  only  the  hydrophones 
actually  present  in  the  HF  array  segment  are  used  in  beamforming,  then 
large  errors  in  the  beam  noise  can  occur  at  tho^e  frequencies  where 
there  are  strong,  spatially  coherent  signals,  such  as  those  due  to 
spectral  lines  in  the  array  tow  ship  noise.  These  errors  have  been 
significantly  reduced  in  the  results  of  this  report  by  using  an 
algorithm  which  forms  estimates  of  the  complex  amplitudes  for  each 
missing  hydrophone  using  the  complex  amplitudes  from  the  four  neigh- 
boring hydrophones.  This  algorithm  is  described  in  Appendix  C and  an 
example  is  presented  which  shows  a reduction  i.n  the  error  in  the  beam 
noise  by  as  much  as  9 dB. 
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APPENDIX  A 

Experiment  and  Data  Processing  Background 

(C)  The  time  series  were  obtained  from  off-line  computations  on 
recorded  hydrophone  data  during  a forty-three  hour  period  from 
1530  hours  on  15  November  to  1030  hours  on  17  November  1974. 

During  this  period,  the  LAMBDA  II  array  system  was  towed  at  3 kts 
by  the  R/V  PACIFIC  APOLLO  along  18  n.m.,  east-west  tracks,  centered 
at  33°  19'  N,  19°  32'  E in  the  Ionian  Basin,  Simultaneously,  a source 
radiating  cw  lines  at  280  Hz  and  40  Hz,  with  source  levels  of  73  and 
155  dB//uPa,  was  towed  at  6 kts.  by  the  R/V  GYRE  along  an  easterly 
track.  The  source  track  and  the  array  tow  site  are  illustrated  in 
the  map  of  Figure  1.1,  ref.  [1].  During  the  observation  period  the 
source  range  decreased  from  148  to  130  n.m.  and  then  increased  to 
174  n.m. 

(S)  The  LAMBDA  II  consists  of  three,  switch-selectable,  array 
segments,  only  two  of  which  are  pertinent  to  this  report.  The  medium 
frequency,  MF , segment  consists  of  64,  equally-spaced,  hydrophone 
groups  yielding  an  array  length  of  812  meters  and  a /2  frequency  of 
60  Hz.  The  high  frequency,  HF , segment  is  152  meters  long  with  a 
2 frequency  of  336  Hz.  It  can  be  viewed  as  consisting  of  67, 
equally-spaced,  hydrophone  groups,  with  three  missing  groups  spaced 
at  quarter-length  intervals  along  the  array.  A detailed  description 
of  the  operati.ng  areas,  the  arrays,  the  on-line  hardware  and  the 
recording  system  can  be  found  in  ref  [1],  along  with  the  results  of 
the  on-line  analysis. 

(C)  The  results  described  i.n  this  report  were  obtained  using  a two- 
dimensional  FFT  beamformer  which  forms  64  beams  and  154  frequency  esti- 
mates at  .2  Hz  intervals  in  each  30  Hz  band,  Hanning  shading  was  used 
in  both  the  time  and  space  domain.  The  3 dB  bandwidth  at  all  frequen- 
cies i.s  .29  Hz.  Values  of  the  3 dB  beamwidths  for  both  array  segments 
at  different  frequencies  and  azimuths  are  presented  in  Table  A.l..... 


MF  Segment 

0.° 

30.° 

60.° 

O 

. 

o 

P- 

30  Hz 

5.2 

6.0 

10.5 

16.5 

60  Hz 

2.6 

3.0 

5.2 

7.7 

HF  Segment 

0.° 

30.° 

60.° 

O 

O 

275  Hz 

3.0 

3.5 

6.0 

9.0 

305  Hz 

2.7 

3.1 

5.4 

8.0 

Table  A.l  (S).  Beamwidths  (3  dB)  (U) 


27 


SECRET 


secret 


(U)  The  individual  samples  of  the  total  beam  power  are  formed  by 
first  averaging  21  consecutive  beamformer  outputs  to  yield  a ten 
minute  average.  An  interpolation  procedure  is  then  used  to  convert 
the  beamformer  variable  to  an  azimuth  variable  after  rotating  in 
azimuth  by  the  amount  required  to  align  the  slightly  different  tow 
track  headings  to  a common  direction. 
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APPENDIX  B 

The  Decomposition  Algorithm  (IJ) 

(U)  In  this  appendix  we  specify  the  algorithm  used  to  decompose 
the  total  beam  power,  B(©,f),  into  the  sum  of  its  spectral  line 
component,  Bs(fl,f),  and  its  background  component,  B^CSjf).  The 
algorithm,  which  is  essentially  the  same  as  that  used  by  Finny 
reference  [3]9lias  four  key  parameters.  These  parameters  are 
identified  and  an  analytical  framework  that  can  be  used  to 
evaluate  their  effect  on  the  decomposition  is  developed. 

(U)  The  decomposition  algorithm  is  based  on  the  principal  that 
the  spectral  lines  can  be  detected  by  comparing  B(fl,f)  against  a 
threshold  that  is  proportional  to  an  estimate  of  the  background 
power.  This  estimate,  which  is  referred  to  here  as  the  threshold 
power  and  is  denoted  by  BT(0,f),  is  computed  by  averaging  B(0,f) 
over  a deleted  neighborhood  of  (0,£),  after  first  discarding  the 
largest  and  the  Kg  smallest  values.  The  purpose  of  discarding 
these  values  is  to  reduce  the  bias  in  B^(0,f)  due  to  the  presence 
of  spectral  lines  that  occur  in  the  neighborhood  of  (0 , f ) . if  a 
line  is  detected  at  <e,f),  the  spectral  component  is  set  equal  to 
the  difference  between  B(0,f)  and  B^(0,f)  and  the  background  component 
is  set  equal  to  the  threshold  power.  If  no  line  is  detected  at  (0,f), 
the  spectral  line  component  equals  zero  and  the  background  power  equals 
the  total  power.  In  either  case,  the  total  power  is  the  sum  of  the 
spectral  component  and  the  background  component  for  all  values  of  (0,f). 

(U)  A precise  specification  is  obtained  as  follows.  Let  B^,  Bg£ 
and  be  the  total  power,  the  spectral  line  power  and  the  back- 

ground power  at  the  ith  frequency-azimuth  resolution  cell 
For  a fixed  resolution  cell,  (©0,£0),  let  RQ  = f (0£,fi);  i=l  , . . . , MJ 
denote  a deleted  neighborhood  of  (0o,fo).  Then  the  threshold  power 
at  ®To’  defined  as 

i=M-KL 

(B.l)  BTo  = ( y~  b'  )/(M-(Kl+Ks)) 


i=Ks+l 

f . . * 

where  the  are  the  order  statistics  0f  B^,  i=l , . . . ,M.  Finally, 

- - ■ — — 

The  B ' are  the  reordering  of  B^  such  that  Bj  £ B2^  ...  ^ bjq 


29 


SECRET 

(This  pagu  is  UNCLASSIFIED) 


SECRET 


tant  c,  the  spectral  line  component 
at  (0o,fQ)  are  defined  by, 

^ Bo  ^ c BTo 
otherwise 


if  B0  ^ c Bt 
otherwise 

(U)  There  are  four  basic  parameters  in  the  decomposition  algo- 
rithm: the  threshold  constant  (c) , the  shape  of  the  region  R0, 

the  number  of  points  (M)  in  Ro  and  the  numbers  of  large  and  small 
values  that  are  discarded,  (Kj)  and  (Kg).  The  spectral  line  and 
the  background  components  described  in  Section  2 were  obtained 
with  c=1.8,  K^  = Kg  = 1,  M=8  and 

(B.3)  R0  = { (6o.fi)  ; i = t 1,  I3.i5.i7  l 

Thus,  the  results  in  this  report  were  obtained  with  B^0  computed 
as  an  average  of  six  values  from  a 2.93  Hz  frequency  band  on 
each  beam  after  first  discarding  the  largest  and  the  smallest 
values . 


(U) for  a given  threshold  cons 
and  background  components 


(B.2a)  Bso=\Bo-BTo 

10 


and 


(B.2b)  Bbo=|BTo 


(U)  In  order  to  rigorously  evaluate  the  effect  of  these  parameters 
on  the  decomposition,  it  is  necessary  to  specify  the  probability 
law  on  the  total  beam  power  at  each  frequency-azimuth  resolution 
cell.  It  is  then  possible,  at  least  in  principal,  to  estimate  the 
number  of  actual  spectral  lines  present  in  B(6  , f)  in  terms  of  the 
number  of  spectral  li.ries  observed  in  Bg(fl  , f) . In  particular,  it 
can  be  shown  that  if  Pg  is  the  fraction  o f the  (©, f)  pi  ane  occupied 
by  actual  spectral  ines  and  Pg  is  the  fraction  of  the  (0,f)  plane 
for  which  B (6,f)  is  non-zero  then 

o 

(B.4)  Ps  - (Pg/PD)  (1-Pf/Ps)/(1-Pf/Pd) 


where 

and 

PF  ' 

are 

! the 

probabi 

.li 

defined 

t>y. 

(B.5a) 

PD 

= 

Prob 

[ Z^c 

i 

(B.5b) 

PF 

= 

Prob 

[ Z>yc 

i 

where 

Z = 

Bo 

/ 

BTo 

ties  of  detection  and  false  alarm 

spectral  line  at  (0o,fo)] 
no  spectral  line  at  (®0,fn)] 
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(U)  Furthermore,  the  fraction  of  the  non-zero  portion  of  Bg(©  ,£) 
that  corresponds  to  actual  spectral  lines  and  the  fraction  that 
results  from  the  statistical  variations  in  B(0,  f)  when  no  lines 
are  present  are  equal  to  Pg/Pg  and  Pp  (l-“^)/Ps  respectively, 
where  Pg  is  given  by  equation  B.4. 

(U)  Idaally,  the  parameters  of  the  decomposition  algorithm 
should  be  chosen  so  that  Pp  is  close  to  unity  and  Pp  is  very 
small.  Then  P will  be  approximately  equal  to  and  Bs(©,f) 
will  be  an  accurate  representation  of  the  spectral  lines  actually 
present  in  B(0,f).  The  effect  of  increasing  the  threshold  con- 
stant is  to  reduce  both  Pp  and  Pp,  and  hence  to  reduce  the  number 
of  lines  observed  in  Bg(9,f).  This  effect  has  been  observed  in 
Figure  3.4.  The  computation  required  to  rigorously  evaluate  this 
effect,  as  well  as  the  effect  of  the  remaining  parameters,  is 
beyond  the  scope  of  this  report.  Insead,  we  first  derive  bounds 
on  the  probabilities  of  detection  and  false  alarm  which  can  be 
used  as  the  basis  for  further  study.  Secondly,  a weak  upper  bound 
on  the  probability  of  false  alarm  is  evaluated  numerically  by 
assuming  that  the  beam  power  is  distributed  as  a chi-squared 
random  variable  when  no  spectral  lines  are  present.  Finally, 
the  effect  of  the  shape  of  Rq  is  illustrated  with  several  examples. 


(U)  The  bounds  on  the  probabilities  of  detection  and  false  alarm 
are  obtained  as  follows.  Let  pn  be  the  probability  that  exactly 

r.nntain 
O 


n resolution  cells  in  R„  contain  power  from  spectral  lines  and 


let  PD 

In  and 

PFln  b 

e the  conditional  probabilities 

and  fa 

Ise  al 

arm  aef 

ined  by, 

(3.6a) 

PDln 

= Prob 

[zz  C 

1 spectral 

line  at  (©( 

V. 

1 spectral 

power  at  n 

(B  ,6b) 

PFln 

= Prob 

fz  >,  c 

■ no  spect: 

ral  line  at 

L 

1 spectral 

power  at  n 

Then  P 

p and 

Pp  can 

be  writt 

en  as 

(B  .7a) 

PD  = 

PD|n  Pn 

(B  .7b) 

PF  = 

E 

PFJ.n  Pn 

Next , 

it  is 

noted  t 

hat  both 

the  condti' 

i n a i.  probabi 

,1? 


and  I 

in  RoJ 

(S  f0)  ancTj 
cells  in  R„  i 


"J 


les  rDln  and 

Pp^n  are  decreasing  functions  of  n,  since  increasing  the  number 
of  spectral  lines  in  R0  can  only  increase  b^,  and  hence,  must 
reduce  the  probability  that  BQ  exceeds  c B.p0.  It  follows  from 
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this  observation  and  equation  (B-7b)  that  Pp  satisfies, 

(B.7c)  PflKL  PKL^  PF  ^ PF1o  PKL  + PFIKL+1(‘1“  PKL^ 
where 


(B  .9) 


Pn 


n=o 


with  a similar  ir.eqality  for  P^. 

(U)  The  bounds  on  Pp  in  the  inequality  (B.8)  should  be  reasonably 
tight  when  Kj  is  large  enough  so  that  PKL  is  close  to  unity.  When 
this  is  the  case,  the  second  term  on  right-hand-side  or  (B.8)  can 
be  neglected.  Furthermore,  Pp|0  and  Ppi^^  should  be  approximately 
equal, since  the  beam  power  in  the  resolution  cells  containing 
spectral  lines  should  be  discarded  before  computing  the  average. 
Thus,  for  large  , Pp  = PpjQ  Ppp  aa^  a similar  argument  can  be 
used  to  conclude  that  = Pp,j0  Ppp- 


(U)  To  numerically  evaluate  the  upper  and  lower  bounds  for  Pp  and 
Pp,  it  is  necessary  to  specify  pn  for  n=o,...,  and  to  compute 
the  conditional  probabilities  for  n=o,  and  Kj+1.  These  tasks 
require  the  specification  of  the  probability  law  on  the  locations 
and  the  power  in  the  spectral  lines  which  is  generally  unknown 
apriori.  Furthermore,  when  K,  and  Kg  are  not  zero,  the  computa- 
tion of  the  conditional  probabilities  is  complicated  by  the  fact 
that  the  probability  law  on  B,p0  is  induced  by  the  probability  law 
on  the  order  statistics  of  the  B^, as  indicated  in  equation  B.l. 

To  alleviate  these  difficulties,  we  first  note  that  since  the. 
conditional  probabilities  are  decreasing  function  of  n,  the  right- 
hand-side  of  B.7c  is  bounded  above  by  Ppi0  and  the  left-hand-side 
is  bounded  below  by  pQ  Pp^0»  Thus,  Pp  also  satisfies 


(B.10) 


PFLo  ^ 


P i 


FLo 


with  a similar  inequality  for  P^.  These  bounds,  although  weaker 
than  those  in  B«7c,  require  only  the  specification  of  po  and  the 
computation  of  the  condtional  probabilities  given  no  spectral 
lines  in  RQ . Finally,  it  is  noted  that  when  and  Kg  are  ap- 
proximately equal  and  no  spectral  lines  occur  in  R , the  proba- 
bility law  on  - b/  Tu’d  be  approximately  the  same  as  that 
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(U) induced  directly  by  the  probability  law  cn  the  B^.  Thus,  for 

equals  Kg  , the  conditional  pro>  ibilities  Pp^Q  and  P^.  can  be 
approximated  by  viewing  B^o  as  an  average  over  M-(Kp+KgJ  points 
and  otherwise  neglecting  the  effe<  " of  discards. 

(U)  The  upper  bound  on  Pp  is  obtained  by  specifying  the 
probability  law  on  the  total  beam  power  given  no  spectral  lines 
in  RQ  and  computing  Pp.  . For  specific  results  we  assume  that 
the  B^  are  proportional  to  independent,  identically  distributed, 
chi-squared  random  variables  with  2r  degrees  of  freedom,  where  r 
is  the  number  of  samples  in  the  time  average  used  to  obtain 
This  assumption  is  consistent  with  a Gaussian  noise  field, 
independent  time  samples  of  beam  power  and  points  in  R suffi- 
ciently separated  so  that  the  B^  are  uncorrelated.  It  follows 
from  this  assumption  that  B^,0  is  proportional  to  a chi-squared 
random  variable  on  2r(M-(Kp+Kg) ) degrees  of  freedom  and  that  Z 
is  F distributed  with  2r  and  2r (M-( K^+Kg ) ) degrees  of  freedom. 
Thus,  Ppi0  can  be  evaluated  for  different  threshold  constants 
from  tables  of  the  F distribution.  In  particular,  for  21  time 
samples  and  the  parameters  used  in  the  decomposition  algorithm 
Pp  is  bounded  above  by  Pp^0  = .004. 

(U)  We  conclude  this  appendix  with  a discussion  of  the  effect 
of  the  shape  of  RQ  on  the  decomposition.  This  effect  can  be 
seen  by  viewing  B^.(0  ,f)  as  a smoothed  version  of  B(0  ,f).  Then 
the  spectral  component  is  that  part  of  B(0  ,f)  that  exceeds 
BT(©  ,f)  by  the  threshold  constant  and  B^( 0 ,f)  is  the  differ- 
ence between  B(Q  ,f)  and  B (0,f).  In  this  context,  the  shape 
of  R()  determines  a two-dimensional  filter  on  the  (ft,f)  plane* 
and  its  effect  can  be  examined  by  considering  the  smoothing 
properties  of  this  filter  relative  to  the  (0,f)  variations  in 
B(0,f). 

(U)  The  plots  of  the  total  beam  power  in  Section  2 suggest 
that  the  variations  in  B(©  ,f)  can  be  visualized  in  terms  of 
three  components:  a component  that  is  smooth  in  both  0 and  fj 
which  represents  the  shipping  noise  that  is  not  resolved** 


* Strictly  speaking  the  associated  filter  is  data  dependent 
due  to  the  Kj +Kg  discards. 

**For  a given  noise  field  the  relative  contributions  of  the 
the  three  components  depends  on  the  array  length  and  the 
analysis  bandwidth. 
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(U)  in  © and  f;  a component  that  is  smooth  in  f but  not  in  ©;, 

which  represents  the  frequency  continuous  portion  of  the  spectra 
from  noise  sources  that  are  resolved  in  © and  finally,  a com- 
ponent that  is  smooth  in  neither©  or  f,  which  represents  the 
spectral  lines  from  noise  sources  that  are  resolved  in  both 
© and  f . 

(U)  The  shape  of  the  region  Rq  used  to  obtain  the  results  in 
the  report,  see  equation  B-3,  represents  one  extreme,  since  this 
shape  corresponds  to  a filter  with  frequency  extent  but  no  azimuth 
extent.  The  effect  of  this  filter  is  to  smooth  the  frequency 
variations  in  B(@  , € 'l  while  preserving  the  azimuthal  variations. 

As  a result,  the  decomposition  algorithm  tends  to  assign  only  the 
spectral  lines  to  Bg(©  ,f). 

(U)  At  the  other  extreme,  consider  a region  which  has  azimuthal 
extent  but  no  frequency  extent.  Such  a filter  would  smooth  the 
azimuthal  variations  in  B(©  ,f)  while  preserving  the  frequency 
variations.  With  this  filter,  the  decomposition  algorithm  tends 
to  assign  both  the  spectral  lines  and  the  frequency-continuous 
components  to  Bg(0  ,£).  This  behavior  is  illustrated  in  Figures 
B-l  and  B-2  which  show  Bg(0  ,f)  and  B^(0,f)  for  sample  one  in  the 
MF  band  with  c=1.8,  M=14,  KL=Kg=2  and 

Ro  “ $ (®i  > ^ • • • > * 7 \ 

(U)  The  comparison  of  these  plots  with  Figures  2.5  and  2.7  illus- 
trates the  difference  in  the  results  for  the  two  extreme  choices  of 
the  shape  of  RQ . 

(U)  As  a final  example,  consider  the  region, 

*0  =i(ek,fi)  ; i=i  1,  1 3,  ± 5,  1 7 ; fe  = -1,0,  1 \ 

(U)  This  region  has  a small  azimuth  extent  and  the  same  frequency 
structure  as  the  region  of  equation  B„3.  As  a result,  the  correspond- 
ing filter  smooths  both  the  6 and  f variations  in  B(®,f).  The  spectral 
component  for  sample  one  in  the  MF  band,  obtained  usi  ng  this  region 
with  c=i.8,  M=24  and  K^=Kg=3>is  shown  in  Figure  B.3.  This  plot  is  essen- 
tially the  same  as  Figure  2.5,  except  that  some  of  the  spectral  lines 
that  are  closely  spaced  in  azimuth  and  are  not  totally  separated  in 
Figure  2,5,  are  totally  separated  in  F'igurc  B.3.  This  example  suggests 
that  some  azimuthal  extent  in  RQ  can  improve  the  resolution  of  the 
spectral  lines  while  keeping  the  frequency-continuous  portion  of  the 
spectra  in  8.^(0.  f) 
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AZIMUTH  (DEGREES) 


Kig.  U.  1 (C)  — Spectral  line  component;  azimuth  smoothing  only, 
56-7-1  dl!//uPa;  ME  hand,  sample  one  (U) 
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AZIMUTH  (DEGREES) 

Fig.  B.2  (C)  — Background  component;  azimuth  smoothing  only, 
5 ( j- 5 7 dB//uI’a;  ME  band,  sample  one  ( IJ ) 
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Fig.  b.3  (C)  — Spectral  line  component;  both  azimuth  and  , equency  smoothing, 
56-74  dB//uPa;  MF  band  sample  one  (U) 
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APPENDIX  C 

An  Algorithm  for  Reducing  the  Effect  of  Missing  Hydrophones  (U) 

(U)  The  HF  array  segment  can  be  viewed  as  an  array  of  67,  equally- 
spaced,  hydrophones*  with  three  missing  hydrophones  located  at  ap- 
proximately quarter  length  intervals  along  the  array.  The  effect 
of  these  missing  hydrophones  is  to  increase  the  peak  sidelobe 
levels  of  the  beam  response  to  a spatially  coherent  signal 
to  within  24  dB  of  the  maximum  beam  value.  This  does  not  significantly 
alter  the  azimuthal  distribution  of  the  beam  power  at  frequencies 
where  no  strong  coherent  signals  are  present.  However,  at  those 
frequencies  with  coherent  signals  that  yield  beam  signal-to- 
noise  ratios  in  excess  of  24  dB,  the  azimuthal  distribution  of 
the  beam  power  will  be  dominated  by  the  sidelobe  levels  over 
the  full  180  degree  azimuthal  sector.  In  this  Appendix,  we 
specify  the  algorithm  that  was  used  to  reduce  the  sidelobe 
levels  due  to  the  missing  hydrophones  and  show  an  example  illus- 
trating the  difference  in  the  beam  response  obtained  using  the 
algorithm  and  that  obtained  by  beamforming  using  only  the  hydro- 
phones that  are  actually  present. 

(U)  The  algorithm  computes  an  estimate  of  the  complex  amplitude 
for  each  missing  hydrophone  using  the  complex  amplitudes  from 
the  four  neighboring  hydrophones.  The  amplitude  estimate  is 
computed  as  an  average  of  the  amplitudes  of  the  four  neighbor- 
ing hydrophones.  The  estimate  of  the  phase  factor  is  computed 
by  first  forming  two  estimates,  each  involving  three  of  the 
four  neighboring  hydrophones,  and  then  combining  these  esti- 
mates to  obtain  the  average  phase  factor.  Finally,  conventional 
frequency-domain  beamforming  is  applied  to  the  augmented  set  of 
complex  amplitudes. 

(ll)  A detailed  specification  of  the  estimation  procedure  is  as 
follows,^  Let  be  the  complex  amplitude  of  the  ic ’ hydrophone 
and  let  Z-  be  the  estimate  of  the  complex  amplitude  of  the 
missing  hydrophone  iQ . Then  the  amplitude  JZ£0\  is  computed  as, 

(c-l)  lZio\  = <|zio_2|  + 1 Zio-jl  + Uio+1l+  |zio+2l  )/4 

and  the  phase  factor,  exp$j©^0^  is  computed  as, 

(C-2a)  eXP{  jGioJ  « <Wm+Wp)/(  |wm  * Wp  | ) 


*The  "hyd rophones"  referred  to  in  this  appendix  are  actually 
groups  of  individual  hydrophone  elements. 
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(U)  where 


(C-2b) 

Wm 

= <zio-2 

Zio  + l 

Zi^-1)/  1 Zio~2 

■'io+1  Zio-l 

(C-2e) 

WP 

^Zio+2 

z* 

^io+1 

Zio-l ^ ! 1 zio  + 2 

1 

o 

• rA 

S3 

*— < 
+ 

o 

• 

S3 

(U)  It 

is 

noted  that 

the 

computations  involve  only  the 

amplitudes  and  their  magnitudes.  This  is  done  to  avoid  the 
2TT  ambiguity  that  arises  in  computations  on  the  phases  asso- 
ciated with  the  complex  amplitudes. 

(U)  The  algorithm  can  be  explained  in  hueristic  terms  as 
follows.  Consider  first,  those  frequencies  where  a strong 
spatially  coherent  signal  is  present..  At  these  frequencies,  the 
amplitude  distribution  in  the  vicinity  of  the  missing  hydrophone 
will  be  approximately  constant  and  the  phase  distribution  will 
be  nearly  linear.  When  this  is  the  case,  the  estimate  will 
approximate  the  true  complex  amplitude  at  the  missing  hydrophone 
and  the  sidelobe  levels  due  to  the  strong  coherent  signal  will 
be  significantly  reduced,  To  see  this,  assume  that  can  be 
represented  as 

(C-4)  Zi  = 1 Z±\  exp  { j(Ki+  * J \ 

where  Ki  is  the  phase  due  to  the  coherent  signal,  with  K depending 
on  both  the  frequency  and  the  incident  angle,  and  ^ is  the 
contribution  due  to  noise.  Substitution  of  (C-4)  into  (C~2),  and 
some  algebraic  manipulation,  yields, 

(C-5a)  exp  j j 0 = exp  { j ( KiQ  + 

where 


(C-5b)  V>c  - (Vi-2  ♦ Vi.p/2  * «■ 

<C-5c)  •»  = ) 0 i£  cos  l C i-2  - P it2>/2  * <Pitl  >/  ° 

otherwise 


• JL.—  JL 

10  il 

ir  c 


For  large  signal-to-noise  ratios,  ot.  will  almost  always  be  zero, 
will  generally  be  small,  and  hence,  exp  { jQ^will 
approximate  the  true  phase  factor.  Furthermore,  the  amplitude 
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(U) of  the  estimate,  which  is  computed  according  to  equation  (C-l), 
will  also  approximate  the  true  amplitude.  Thus,  the  beam  response 
obtained  by  beamforming  on  the  augmented  set  of  complex  amplitudes 
will  have  sidelobe  levels  that  approximate  those  that  would  result 
if  the  missing  hydrophones  were  present. 

(U)  Next  consider  those  frequencies  where  no  strong  spatially 
coherent  signal  is  present.  At  these  frequencies,  both  the 
amplitude  and  the  phase  factor  of  the  estimate  may  differ  signif- 
cantly  from  the  true  values  at  any  given  time.  The  amplitude  of 
the  estimate  will  be  roughly  equal  to  the  mean  amplitude  due  to 
the  noise  alone  and  the  phase  of  the  estimate,  relative  to  the 
true  value,  will  be  approximately  distributed  as  a uniform 
random  variable.  However,  the  correlation  between  the  phase 
factor  of  the  estimate  and  those  corresponding  to  the  neighbor- 
ing hydrophones  will  be  small.  Thus,  the  beam  response  obtained 
using  the  augmented  set  of  compelx  amplitudes  will  again  be 
approximately  the  same,  in  the  average,  as  is  if  the  missing  hydro- 
phones were  present, 

(U)  A rigorous  performance  evaluation  of  the  algorithm  requires 
the  specification  of  the  probability  laws  on  both  the  signal  and 
the  noise  fields  and  the  computation  of  the  mean  sidelobe  level 
obtained  using  the  estimation  procedure.  This  computation,  wnich 
is  complicated  by  the  fact  that  the  estimation  procedure  is  non- 
linear, is  beyond  the  scope  of  this  appendix.  In  lieu  of  this 
evaluation,  we  illustrate  both  the  need  for  and  the  effect  of  the 
algorithm,  by  comparing  the  beam  power  from  sample  one  on  the  HF 
array  segment  obtained  using  the  algorithm,  against  that  obtained 
by  beamforming  using  only  the  hydrophones  that  are  actually 
present.  The  plots  of  the  total  beam  power  for  both  cases 
are  illustrated  in  Figures  C-l  and  C-2,  where  the  scale  on  beam 
power  in  each  plot  ranges  from  .38  dB  down  from  the  peak  value 
of  the  280  Hz  signal  line  to  20  dB  down  from  the  peak^,  in  3 dB 
steps  per  shade. 

(U)  The  plot  of  the  total  beam  power  obtained  without  using 
the  algorithm,  Figure  C-l,  illustrates  the  effect  of  the  missing 
hydrophones  on  the  total  beam  power.*  At  several  irequencies,  the 
sidelobe  levels  associated  with  strong  coherent  signals  completely 
dominate  the  beam  power  over  the  full  180  degrees.  This  is  most 
evident  at  the  280  Hz  signal  frequency,  although  there  are  other 
frequencies  where  the  array  tow  ship  noise  is  sufficiently  large 
to  mask  the  azimuthal  distribution  of  the  beam  power  due  to  noise 
sources  other  than  the  array  tow  ship. 

*The  beam  power  in  Figure  C.l,  is  obtained  using  Hanning  shading 
on  64  hydrophones  with  zero  coefficients  at  the  positions  of  the 
missing  hydrophones. 
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(U)  Figure  (C-2)  illustrates  the  beam  power  obtained  by 
including  the  estimates  of  the  complex  amplitudes  in  the  input 
to  the  beamformer.*  A comparison  of  Figures  (C-2)  and  (C-l), 
shows  that  the  sidelobe  levels  associated  with  the  strong  co- 
herent signals  have  been  reduced  significantly.  This  is 
particularly  evident  in  the  low  noise  directions  (in  the  vicinity 
of  £ 50  degrees),  where  the  beam  noise  is  essentially  uncontami- 
nated by  the  side  lobes  from  the  280  Hz  signal  and  the  array  tow 
ship  noise, 

(U)  A more  detailed  comparison  of  the  effect  of  the  algorithm  can 
be  seen  in  Figures  (C-3)  and  (C-4)  where  we  have  plotted  the  azimu- 
thal distribution  of  the  beam  power  at  280  Hz  and  at  276.6  Hz. 

The  continuous  line  curve  in  each  plot  represents  the  beam  power 
obtained  by  neglecting  the  missing  hydrophones  and  the  dotted 
curve  represents  the  beam  power  obtained  using  the  algorithm. 

The  beam  power  in  both  plots  is  normalized  to  the  maximum  beam 
power  at  the  280  Hz  signal  line. 

(U)  The  azimuthaL  power  distributions  at  ,280  Hz,  illustrated  in 
Figure  C-3,  shows  that  the  effect  of  the  algorithm  is  to  signifi- 
cantly reduce  the  beam  power  at  all  azimuths  except  those  where 
spatially  coherent  signals  are  actually  present.**  This  reduction 
is  most  dramatic  in  the  low  noise  directions  where  the  two  curves 
differ  by  as  much  as  9dB, 

(U)  The  azimuthal  power  distributions  at  276,6  Hz,  Figure  C-4, 
shows  that  even  without  the  280  Hz  signal  line,  the  reduction  in 
the  beam  noise  at  azimuths  where  no  dominant  coherent  source  is 
present  can  be  significant.  The  difference  between  the  two  curves 
is  somewhat  less  than  that  in  Figure  C-3,  except  near  -50  degrees 
where  the  two  differences  are  comparable. 


* Three  hydrophones  were  dropped  from  the  ends  of  the  augmented 
set  so  that  an  FFT  algorithm  on  64  points  could  be  used  for 
beamforming. 

**Less  than  .6  dB  of  the  difference  between  the  two  curves  can 
be  attributed  to  the  fact  that  the  scale  for  each  curve  is  rela- 
tive to  the  peak  of  the  280  Hz  signal  line. 
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Fig.  C.3  (C)  — Azimuthal  power  distribution  at  275  Hz  (U) 
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UNITED  STATES  GOVERNMENT 

Memorandum 

7100-034 

DATE:  26  February  2004 

REPLY  TO 

ATTN  OF:  Burton  G.  Hurdle  (Code  7103) 

SUBJECT:  REVIEW  OF  REF  (A)  FOR  DECLASSIFICATION 

TO:  Code  1221.1 

REF:  (a)  “Preliminary  Results  of  an  Analysis  of  Beam  Noise  in  the  Mediterranean”  (U), 

R.  Heitmeyer  and  G.  Long,  Acoustics  Division,  NRL  Memo  Report  3832, 
September  1978  (S) 


1 .  Reference  (a)  is  a report  of  the  preliminary  results  of  an  experiment  to  measure  beam  noise  in 
two  frequency  bands  (30-60  and  275-305).  The  results  indicate  that  the  noise  is  primarily 
due  to  shipping 


2.  The  technology  and  equipment  of  reference  (a)  have  long  been  superseded.  The  current 
value  of  these  papers  is  historical 

3.  Based  on  the  above,  it  is  recommended  that  reference  (a)  be  declassified  and  released  with  no 
restrictions 


CONCUR: 

3 / ) ) 

E.R.  Franchi  Date 

Superintendent,  Acoustics  Division 


CONCUR: 

Tina  Smallwood  Date 

NRL  Code  1221.1 


DEPARTMENT  OF  THE  NAVY 

OFFICE  OF  NAVAL  RESEARCH 
875  NORTH  RANDOLPH  STREET 
SUITE  1425 

ARLINGTON  VA  22203-1995 

IN  REPLY  REFER  TO: 


5510/1 

Ser  32 1OA/0 11/06 
31  Jan  06 


MExMORANDUM  FOR  DISTRIBUTION  LIST 

Subj : DECLASSIFICATION  OF  LONG  RANGE  ACOUSTIC  PROPAGATION  PROJECT 
(LRAPP)  DOCUMENTS 

Ref:  (a)  SECNAVINST  5510.36 

Enel:  (1)  List  of  DECLASSIFIED  LRAPP  Documents 

1 . In  accordance  with  reference  (a),  a declassification  review  has  been  conducted  on  a 
number  of  classified  LRAPP  documents. 

2.  The  LRAPP  documents  listed  in  enclosure  (1)  have  been  downgraded  to 
UNCLASSIFIED  and  have  been  approved  for  public  release.  These  documents  should 
be  remarked  as  follows: 

Classification  changed  to  UNCLASSIFIED  by  authority  of  the  Chief  of  Naval 
Operations  (N772)  letter  N772A/6U875630,  20  January  2006. 

DISTRIBUTION  STATEMENT  A:  Approved  for  Public  Release;  Distribution  is 
unlimited. 

3.  Questions  may  be  directed  to  the  undersigned  on  (703)  696-4619,  DSN  426-4619. 


^ ^ 
BRIAN  LINK 
By  direction 


Subj : DECLASSIFICATION  OF  LONG  RANGE  ACOUSTIC  PROPAGATION  PROJECT 
(LRAPP)  DOCUMENTS 

DISTRIBUTION  LIST: 

NAVOCEANO  (Code  N121LC  - Jaime  Ratliff) 

NRL  Washington  (Code  5596.3  - Mary  Templeman) 

PEO  LMW  Det  San  Diego  (PMS  181) 

DTIC-OCQ  (Larry  Downing) 

ARL,  U of  Texas 

Blue  Sea  Corporation  (Dr. Roy  Gaul) 

ONR  32B  (CAPT  Paul  Stewart) 

ONR  3210A  (Dr.  Ellen  Livingston) 

APL,  U of  Washington 

APL,  Johns  Hopkins  University 

ARL,  Penn  State  University 

MPL  of  Scripps  Institution  of  Oceanography 

WHOI 

NAVSEA 

NAVAIR 

NUWC 

SAIC 
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